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Abstract 

Here we revisit the bag phenomenology of the deconfining phase transition to replenish 
it by introducing systematically the bag surface tension. Comparing the free energy of such 
bags and that one of the strings confining the static quark-antiquark pair, we express the 
string tension in terms of the bag surface tension and thermal pressure in order to estimate 
the bag characteristics using the lattice QCD data. Our analysis of the bag entropy den- 
sity demonstrates that the surface tension coefficient is amazingly negative at the cross-over 
(continuous transition). This approach allows us to naturally account for an appearance of a 
very pronounced maximum (observed in the lattice QCD simulations) of the entropy of the 
bound static quark-antiquark pair. The vicinity of the (tri)critical endpoint is also analyzed 
to clarify the meaning of vanishing surface tension coefficient. 
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1 Introduction 

One of the key physical quantities provided by the lattice quantum chromodynamics (LQCD) 
is the free energy of static quark-antiquark pair F qq -(T, L) as a function of the temperature 
T and separation distance L being extracted from the Polyakov line correlation in a colour 
singlet channel. The linear L-dependence of F qq (T, L) discovered at large L and low temper- 
atures naturally explains the colour confinement. On the other hand at higher temperatures 
the linear L-dependence of F qq (T, L) disappears signaling the Debye screening and an advent 
of deconfined phase. This phenomenologically transparent picture resembles the confining 
string model [1] which is fully adopted by the LQCD community [2, 3, 4]. 

Here we develop another general view of the confinement phenomenon dealing entirely 
with an idea of quark-gluon (QG) bag with nonzero surface tension. In this approach the 
importance of the surface tension concept was realized long ago [5, 6], but only recently the 
surface tension of large QG bags was consistently included into the statistical description of 
the QG plasma equation of state [7, 8, 9, 10, 11]. It turns out quite suitable to formulate the 
analytically solvable statistical models for the QCD tricritical [7] and critical [11] endpoint, 
and to push forward an idea of the finite width model of QG bags [8, 9, 10]. Clearly such 
a development improves our understanding of the QG plasma equation of state and brings 
it forward to a qualitatively new level of realism by establishing the Regge trajectories of 
heavy/large bags both in a vacuum and in a medium [8, 9, 10] using the LQCD data. 
However, to make the approach quantitatively informative one needs to establish the value 
of the surface tension in the whole (maximally possible) range of temperature and baryonic 
chemical potential. Unfortunately, as for now the LQCD cannot provide us with such an 
information. To resolve this problem we are going here to ascertain the phenomenological 



relation between the string tension of a tube confining the static quark-antiquark pair and 
the surface tension of the QG bags in order to study the bag thermodynamics. We believe the 
concept of bag surface tension looks more adequate just at high temperatures. Another major 
(and closely related) task of this study is an investigation of the (tri)critical endpoint vicinity 
of the QCD phase diagram to clarify the meaning of vanishing surface tension coefficient. 



2 Free Energy of Elongated Cylindrical Bag 

The free energy of large almost spherical QG bag can be cast as [7, 8, 9, 10, 11] 

F sp (T,V) = - Pv (T)V + 2a surf (T)S(V) + Tr In [^] . (1) 
Here p v (T) is the thermal pressure inside a bag, (J surf 

(T) is the temperature dependent 

2 

surface tension coefficient, S(V) ~ V» is the mean surface of the QG bag, while the last 
term on the right hand side of (1) is the Fisher topological term [12] which is proportional 
to the Fisher exponent r = const > 1 [7, 11]. Vq is a normalization constant with the 
dimension of volume. Note that this is the standard parameterization of the free energy 
of large physical clusters which is successfully used in the Fisher droplet model [12], in 
the statistical multifragmentation model [13, 14] and in describing the free energy of large 
geometrical clusters of the 2- and 3-dimensional Ising model [15, 16] and of percolation 
clusters [17]. Such a free energy parameterization turns out very efficient in studying the 
critical point of realistic gases [12, 18]. It was applied to many different systems with the 
different extents of success including a nuclear multifragmentation both in infinite [13, 14, 
19, 20, 21] and in finite systems [21], a nucleation of real fluids [22] and the compressibility 
factor of real fluids [23]. 

In principle, besides the bulk (~ V) and surface parts (~ S{V)) the free energy (1) could 
include the curvature part as well, which may be important for small hadronic bubbles [6] 
or for cosmological phase transition study [24]. We stress, however, that as usual the critical 
properties of the statistical models are defined by the infinite bag, therefore, including a 
curvature term of any sign in (1) could affect the thermodynamic quantities of such models at 
(tri)critical endpoint only [7, 11] (see also below). If the curvature term is of a real importance 
for the cluster models discussed here, then it should also show itself at the (tri)critical points 
of many systems described by free energy of Eq.(l) [12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23], 
but this is not the case (see an extended discussion in [7]). Keeping in mind this argument 
we omit the curvature part of bag free energy as well. 

Using the thermodynamic identity 

Ps P = - (££) r = Pv (T) -2. SW/ (T)« - TL , (2) 

one finds the pressure of the spherical bag from (1). To calculate the free energy of bags of 
more complicated shapes one has to change the second term on the right hand side of Eq. (2) 
to the general Laplace form of surface pressure 



Pgen = Pv(T) - a sur f(T) 



1 1 

R\ R2 



f, (3) 



where R\ and R2 are the main curvature radii of the shape defined locally. Basing on Eq.(3) 
one can find the free energy of an arbitrary shaped bag. Then for the large elongated cylinder 
of the radius R = R\ and the height L 3> R one gets (R2 = 00) 

F cyl (T, L,R) = -J dVp gen = - Pv (T)ttR 2 L + a surf (T)27rRL + Tr In V + f(T) . (4) 



Here the volume independent function f(T) is the integration constant. Comparing Eq.(4) 
with Eq.(l) and applying the same steps to obtain the sphere free energy from Eq.(2), we 
conclude that f(T) = — Tt In Vq for a cylinder just like for a sphere. 

Assuming now that the free energy of cylindrical bag (4) of the radius R and the length 
L 3> R equals to the free energy of colour string of the same size F s t r ~ o s t r L that binds the 
static quark-antiquark pair one can find the desired relation at vanishing baryonic densities: 



a stT {T) = a surf {T) 2ttR - p v (T)irR 2 + ^ In 



ttR 2 L 

Vn 



(5) 



In doing so we match an ensemble of all string shapes of fixed L to a mean elongated cylinder, 
which according to the original Fisher idea [12, 18] and the reliable estimates of the Hills 
and Dales Model [27] represents a sum of all surface deformations of a given bag. Choosing 
sufficiently large radius R and very large height L 3> R (thermodynamic limit) one can see 
that the last term in Eq.(5) vanishes. Clearly, the corrections coming from the Coulomb part 
in the string free energy F s t r or from different parameterization of the Fisher topological term 
in Eq.(l) and in Eq.(5) should also vanish in this limit. Note that the term proportional to 
In L is present in the string free energy F str as well, but it was not analyzed in [1] due to the 
same reason. 

The last result shows that for the thin strings R — > (compared to L) or for the van- 
ishing thermal pressure p v — > there exists a simple interrelation between the colour string 
tension and the surface tension of the QG bag a s t r {T) a sur f(T) 2ttR, but in general Eq.(5) 
determines the temperature dependence of string radius 



R ± 



rf(T) ± ^l rf (T) ~ 



Pv(T) 



(6) 



if the temperature dependences of a s t r {T), a sur f(T) and p v (T) are known. On the other 
hand it is also possible to determine the T-dependence of the surface tension of bags 

<WT) = ^P + \p v {T)R, (7) 

if R(T), a s tr{T) and p v (T) are known. In fact, Eq.(6) already gives us the following radius 
independent inequality for bag surface tension 

a surf (Tf > ^T)«,r(T) , (8) 

7T 

which demonstrates that for a confining string, i.e. for a s t r (T) > 0, the bag surface tension 
can vanish only and only for p v (T) < 0, i.e. only for negative or zero values of thermal 
pressure! The latter is clearly seen, if one substitutes cr sur f(T) — >■ into Eq.(6) and considers 
a str {T) > 0, i.e. 



R- 



Vstr(T) (9) 



ir/ (THo V ~*Pv(T) ' 

which is real for p v (T) < only. 

Eq.(7) allows us to estimate roughly the surface tension at T = 0. Taking the typical 
value of the bag model pressure as p v (T = 0) = — (0.25) 4 GeV 4 , R = 0.5 fm and a str (T = 
0) = (0.42) 2 GeV 2 [25], one finds from (7) that a surf {T = 0) = (0.2229 GeV) 3 + 0.5^ R w 
(0.183 GeV) 3 157.4 MeV fm -2 . Our estimate is larger than the ones which are familiar to 
the astrophysics community [24] but we would like to emphasize it is based on conservative 
parameter values in Eq.(7). Optimizing the radius R and bag constant value we could get 
perfectly suitable magnitude of the surface tension at vanishing temperature. 



3 Thermodynamics of Cylindrical QG Bag 



The above results allow us to tune the interrelation with the colour string model and to study 
the bag surface tension near the cross-over to QG plasma phase. The LQCD data indicate 
that at large R the string tension behaves as [26, 4] 



LQCD 
a str 



R 2 



C. 



(10) 



where Lq and C are some positive constants. Such a behavior in a confined phase can be 
easily understood within the confining string model [1]. Indeed, at low T the string energy 
is proportional to the separation L. As the temperature increases, the flux tube starts to 
oscillate and its length exceeds the separation distance L. Eventually the flux tube travels 
all over the available space before ending on the colour charges. 

In fact, a very similar explanation emerges from the view point of elongated cylinder 
despite the different air of Eq.(10) and Eq.(5). Noting that the weak L dependence of 
the LQCD string tension (10) still has to be accurately explored for the large values of 
quark-antiquark pair separation L, here we, however, would like to study the collapsing 
string tension at fixed L 3> Lq (or better to say, the string melting) and the temperature 
approaching the cross-over temperature T tr from below, i.e. for T — > T tr — 0, which according 
to the LQCD [28] and to the flux tube model [1] are equal. Similar to [1] we assume that in 
the infinite available volume a s t r (T) — > +0 => R — > oo in such a way that (ujk = const ~ 
ln(L/L )) 



a str (T) R K -^oo k > 0, 



(11) 



thereby extending a range for the power k > to study more general case, since the formal 
expressions for thermodynamic functions are valid both for positive and negative values of 
k. The value of constant Uk > is not of crucial importance for us here because we are 
interested in the qualitative analysis while it becomes quite essential for the quantitative 
estimates. 

The effect of outer baryonic charge in the system of quark-antiquark pairs (with zero net 
baryonic charge) can be accounted for by the dependence of string melting temperature Tt r 
on the baryonic chemical potential fj, b - Apparently it leads to the temperature Tj r decreasing 
as a function of ^ as expected by the QCD phenomenology [29]. Thus, in what follows we 
also assume that Eq.(ll) is valid for the non-zero values of Performing our analysis at 
fixed values of we have no need to introduce the particular dependence of all quantities 
(including oj^ in Eq.(ll)) on it. We have to keep in mind only that the temperature T tr may 
change with u b . 

The surface tension coefficient can be affected by ^ > as well, but the main result of 
Eq.(5) remains obviously valid. Now neglecting the last term on the right hand side of Eq.(5) 
for L ^> R one can calculate the thermal pressure (treating Eq.(ll)) as 



(rp\ _ r, O-surf(T) CFstr (T) 
p v {l ) - I £ 



2a 



surf 



<£k 



fc+1 
k 



(12) 



Similarly, neglecting the Fisher topological term in Eq.(3) at L » R (for a cylindrical bag of 
radius R and height L 2> R) one obtains the total bag pressure 



Ptot=Pv(T)-- 



(T) 



R 



7 surf 



(T) 



R 



Cstr 



a 



surf 



a str 



■ (13) 



and its total entropy density 



Stot 



d pt 



ot 



d T 



<^k 



k d CFstr 

d T 



rf + 



1 

k d 



surf 



d T 



k+2 

TT k 



^k 



k d CFstr 

d T 



(14) 



The mechanical stability of the cylindrical bag means an equality of the total bag pressure 
Eq.(13) to the outer pressure, but the thermodynamic stability requires positive value for 
the entropy density (14). To quantify the latter we adopt the following parameterization of 
the string tension for T — > T tr — 



a str (T) = a° str t\ where t = ^—^^+0, (15) 



'-tr 



with a° str > and v > (for example, in the simplest case v = 1 at \i = [1], but the other 
values can also be valid). Then the total entropy density of the cylindrical bag becomes 



Stot 





1 


J V u surf 






\ kT tr t 



i 9 a aur f (k+2)u 
dT nk 



,t v 



Ttr 



(16) 



The model of quark gluon bag with surface tension [7, 11] predicts that everywhere at the 
cross-over line, except for the (tri)critical endpoint, the surface tension coefficient cr sur f is 
non-zero and its derivative Q S !f f is finite at t — > +0. Remembering this requirement one 
finds from Eq.(16) that its first term of the right hand side dominates and, hence, we receive 



Stot -> - 



a str 



i 



, f {T tr ) t^ 1 > , (17) 



kTi 



tr 



which means that at T — > Tf r — the surface tension coefficient must be negative & surf {Ttr) < 

0. Actually, this result brings nothing surprising since the calculations of surface partitions 
for physical clusters [27] and the model of quark gluon bag with surface tension with tricritical 
[7] and critical endpoints [11] predict also that at low baryonic densities the deconfining phase 
transition degenerates in a cross-over just because the surface tension coefficient of large bags 
becomes negative in this region (for more details see next section). Eq.(17) clearly shows 
that the colour string model shares the possibility of negative values of bag surface tension 
coefficient available in the cross-over region. 

However, considering the above results in the context of the LQCD data we should men- 
tion it is unlikely that the current calculations on the finite lattices allow us to claim an 
immutable validity of Eq.(ll). Indeed, an analysis of the quarkonium spectra in the de- 
confined phase teaches they are surviving (not melted) far behind the critical temperature 
signaling that the string tension does not completely vanish at high temperatures [28]. Ap- 
parently, it implies the modification of Eq.(ll) and new phenomenological inputs inevitable. 
The LQCD data for the free energy (and entropy) of the colour string [28] demonstrate an 
extremely fast increase at approaching the cross-over temperature from below although the 
mean free energy of the string (F s t r ) is finite at L — > oo. In order to make these LQCD 
results instrumental and not to modify Eq.(ll) at the scale of finite lattice of the spatial size 
Riat the mean free energy of the string at L — > oo can be subdivided in two parts. One part 
with the probability W(L) corresponds to the strings of infinite values of Fj£ r = a s t r L in the 
limit L — > oo and another part for the strings having the finite values of free energy -F s ° tr , 

1. e. {F str ) = F£t r W + Fg tr {l — W). Then it seems reasonable to assume that the entropy of 
the strings of finite free energy is finite and smooth function of t — > +0. Now to explain the 
behavior of the lattice free energy and entropy of the colour string basing on Eqs.(ll) and 
(15) we have to expect W{L) behaving as W ~ [L ln(L/Lo)] _1 at L — > oo and fixed R. If 
it is well argued to use the above results for R < Ri at and L — > oo to compare them with 
the LQCD data, then using the thermodynamic identity = —S^ r , Eq.(5) and the first 
equality in Eq.(13) one obtains the lattice entropy S^ r at fixed L as 



S ^tr = ~ _ . ^— -T5 ^——W L - (18) 



-irR - stotirR 2 



Q- sur fTrR _ 2p t otTrR 2 

k 0~ str k 0~str 



The colour tube radius for vanishing a str can be found from Eq.(6) as R -> ffsari. It 
helps to simplify the denominator in Eq.(18) and then neglecting the term with the derivative 
of surface tension in the numerator of Eq.(18) (the term with stot is larger) one finally receives 

c oo _^ stotka s t r R T s tot ku k 

u surf u surf 1L 

where at the last step of deriving (19) we used Eq.(ll). 

It is easy to verify (using Eqs.(15) and (17)) that despite the approximations done Eq.(19) 
is nothing more than (-^L). Now it becomes clear that the contribution of this term 
into the total entropy is finite S^ r W < oo for L — > oo and R < R[ a t- Again we see that it 
should be cr sur f(T tr ) < to provide the positive values of S^ r and s to t- 

Besides, it is clear from (19) that with t decreasing to a minimal value and R approaching 
Riat one has S™ r W ~ t v ~ l , i.e. the entropy part of the strings of infinite free energy increases 
quickly for v < 1. A formal extension of this result to R — > oo leads to a divergency of 
S%t r W ~ t u ~ l at t = 0. However, to study the behavior of the lattice entropy with more 
accuracy it is necessary to assume a certain behavior of the string tension in the vicinity of 
t = 0. 

In order to demonstrate the possibility for entropy density Eq.(14) to be divergent at 
t = even for the lattice of finite size we consider the simplest modification of Eq.(15) which 
clarifies the fact of string tension survival at t = 0. It has still a finite magnitude but is going 
down for T > T tr [28] . Actually, such a modification accounts for that the radius of the colour 
string R in (11) can not exceed the lattice size. We choose the simplest parameterization 

^"U + ^-Hr, T^T tr + 0, (20) 

where t = (T tr — T)/T tr , the string tension coefficient at the cross-over temperature a% r > 
is small but non-zero, a ± and v ± > are some non-zero constants. Now using Eq.(20) we 
find the string tension derivative 

dastr = J_i -a-u-t^- 1 <0, T^T tr -0, 

dT T tr \ a + v + {-t) v+ - 1 <0, T^T tr + 0. { ' 

Since the LQCD data show that the string tension coefficient is monotonically decreasing 
with T on both sides of T = T tr , whereas the entropy increase (decrease) is very fast for 
T->T tr -0(T-> T tr + 0) [28], we conclude from Eq.(14) and Eq.(21) that such a behavior 
can be provided by — 9 q s ^ L W(L) if the following inequalities are valid o~ > 0, a + < 0, 
<J surf {Ttr) < and ^ ± < 1. In proving this statement one has to account the finiteness of 
& surf (Ttr) and its derivative and the smallness of <7*j r compared to \& sur f 

Thus, the presence of a non-zero surface tension in the entropy density (14), (16) and 
(17) naturally explains also the origin of extremely fast entropy increase close to the cross- 
over temperature for any lattice sizes. This effect was called 'mysterious' in Ref.[29] but we 
note here that a similar behavior should be inherent in the energy density as well since the 
pressure Eq.(3) is finite. In principle, one could choose more sophisticated T-dependence of 
the string tension a s tr in order to describe the finite maximum of the lattice entropy but 
such a task requires, first of all, the LQCD data of very high quality to be conclusive. 

Furthermore, the equations (14), (16)-(21) derived together with Eq.(ll) provide us now 
with a principal possibility to determine the QG bag surface tension along the cross-over line 
in the (p,b — T)- plane directly from the LQCD data. Such results are of vital importance for 
quantitative estimates while dealing with the model of quark gluon bags with surface tension 
[7, 11] and its generalizations [8, 9, 10]. 



4 The (Tri) critical Endpoint Vicinity 



Another special and interesting possibility to be considered is if the surface tension coefficient 
in Eq.(16) vanishes simultaneously with the string tension as 



rf(T) 



a surf ■ tC 



surf 



T^T tr -0, 
T —■ T tr + , 



(22) 



or, in other words, the surface tension vanishes at the cross-over line. As shown in [7, 11] 
such a situation takes place at the tricritical or critical endpoint only, whereas in the example 
of Eq.(17) the surface tension coefficient is considered at the values of baryonic chemical 
potential which are smaller than that of the (tri)critical endpoint. Here we use the surface 
tension in the form introduced in [11] intending to apply it to the analysis of (tri)critical 
endpoint. 

Choosing this simple surface tension parameterization, Eq.(22), we follow the original 
Fisher idea [12] to explain the temperature dependence of surface free energy for = 1. 
According to that idea the surface free energy consists of two terms, the surface energy of 



a bag of volume V as crf ur j 



1 + — 1 2 

V» and the term —Taf ur jT tr Vs which comes from the surface 



entropy o'furfT^Vs in the original model [12]. Note that the surface entropy of a bag of 
volume V counts its degeneracy factor or the number of ways to have such a bag with all 
possible surfaces. This result was generalized to > 1 using the exact solutions of the 
Hills and Dales Model [27] while we consider the interval < < 1. The Hills and Dales 
Model naturally explains the negative values of the surface free energy. They are originated 
by the dominance of bags of complicated (non-spherical) shapes whose number of states 
(exponential of entropy) simply exceeds their suppression by the Boltzmann factor. 

Substituting (22) into (16) and using (15) we get for T — > T tr — in thermodynamic limit 



Stot 



[c- 



surf 



Ttr 



(k+2)v 

TT k 



T str L 



1 

k & 



4.1/- 

str \ 



Ttr 



(23) 



In contrast to (17) the surface tension term in (23) dominates for ( < ^-^v only and it 
implies a~ _ . < 0. Let us note the latter contradicts to critical endpoint existence conditions 

sur j 

found in [11] but such an inequality does not affect the tricritical endpoint existence [7, 8]. 
The entropy density would diverge in this case for £~ < 1 — j£ (which leads to v < k) whereas 
for 1 — | < (~ < ^-j^v (which is true for v > -r^—) the entropy density would vanish at 
the tricritical endpoint. If, however, £~ > then the last term on the right hand side 

of Eq.(23) dominates and, hence, the entropy density is obviously positive for k > and 
v > 0. Under this condition for £ - both the critical and tricritical endpoints may exist and 
the entropy density diverges at t = for v < 

Considering now Eq.(20) as the finite size analog of (15) we have for the entropy density 
in the vicinity of =pi — > +0 



Stot 



tr 
'sir 



k (±C 



Ttr 



sur j 



(fc+2) 

7T k 



Uk 



+ 



surf 
Ttr 



C ± (Tt f 



-l 



(24) 



This equation has more complicated structure being compared to Eq.(23). For > and 
> 0, however, the term which results from the surface tension in (24) (proportional to 
(-ft) u +< ' X ) can not dominate anymore. Nevertheless, the term originated by its derivative 
(the last term in (24)) is dominating for v ± > ^ implying ±af ur jr > 0, since according to 

(21) ±a ± < 0. And as before this possibility, ±crf ur j > 0, is inconsistent with the critical 



endpoint existence [11]. For < the second term on the right hand side of Eq.(24) 
dominates and it leads to the positive value of entropy density (24). 

Clearly, the entropy density Eq.(24) would diverge at t = for either < 1 or < 1 
being essentially different from what we found out at analyzing (23). The important message 
of this comparison is that one has to be very careful at operating with the LQCD results to 
extract a physical information above a cross-over on the phase diagram. 

5 Conclusions 

Here we develop the novel approach to the bag phenomenology of deconfinement using quark 
gluon bag model with surface tension [7, 8, 9, 10, 11] as an effective tool for exploring the 
particular region of phase diagram of strongly interacting matter. Analyzing the free energies 
of cylindrical bags and confining colour strings we find out the quite general alignment 
connecting the string tension with bag surface tension and thermal pressure. This relation 
makes it possible, in principle, to estimate the bag surface tension directly from the LQCD 
measurements. We use the relation derived to study (under the plausible assumptions) the 
string entropy behavior in the cross-over vicinity, i.e. for (almost) vanishing string tension, 
taking into account the finite lattice size in the LQCD simulations. We find out the divergent 
behavior of the entropy density for the finite values of bag surface tension at the cross-over. 
The relation (obtained in Eq.(19) between the cylinder bag entropy density and the lattice 
entropy S^f r W(L) allows us to naturally explain the 'mysterious' [29] maximum of S^f r W(L) 
measured in the LQCD simulations at the cross-over. Moreover, drawing the corresponding 
LQCD results we determine that the surface tension coefficient is amazingly negative at the 
cross-over in conformity with the prediction of quark gluon bag model with surface tension 
[7, 8, 9, 10, 11]. This model says that the only physical reason for transforming the first 
order deconfining phase transition into the cross-over at low baryonic chemical potentials 
is just the negative values of the surface tension coefficient of bags in this region of phase 
diagram. Apparently, the LQCD data give us the quite serious confirmation of this finding. 
Comparing the behavior of entropy density in the vicinity of (tri)critical endpoint for infinite 
and finite systems we find out the different quantitative results although their qualitative 
behaviors are similar. Besides, the estimates obtained could be much more conclusive if rely 
upon the reliable LQCD data at non-zero baryonic densities. Unfortunately, we are still only 
at the starting point on this road and due to this fact the bag model with surface tension 
looks more suitable for phenomenological analysis, at least, today. 

The formulae derived allow us to estimate the bag surface tension at zero temperature 
drawing the standard values of thermal pressure as given by the equation of state of bag 
model. Surprisingly it occurs rather large a sur f(T = 0) ps (0.183 GeV) 3 ps 157.4 MeV 
fm -2 if compared to the estimates received in the other phenomenological models. Surely, 
more accurate LQCD analysis of the bag surface tension will be done in a time to come 
but today we believe our estimate based on the conservative values of parameters provides, 
in a sense, its lower bound which could be practical for studying the early Universe. In 
one of those scenarios [30] the quark matter may survive the 'boiling' period of the early 
Universe if a sur f > (0.178 GeV) 3 . However, the idea of quark matter survival to the present 
days discussed in [31] requires a revision in view of the negative values of surface tension 
coefficient at the cross-over. Moreover, it definitely needs an additional phenomenological 
input since the large QG bags have very short life-time as shown within the models like 
[8, 9, 10] or within the Hagedorn-Mott resonance gas model [32]. 
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